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Abstract 

We provide formulas and algorithms for computing the excess numbers of certain ideals. 
The solution for monomial ideals is given by the mixed volumes of certain polytopes. These 
results enable us to design specific homotopies for numerical algebraic geometry. 

1 Introduction 

Consider a homogeneous ideal I C C [xq, . . . , Xn] generated by the forms Bi, . . . ,Bk of degrees 
Pi, . . . ,Pk. Let /i, . . . , be forms of X of degrees di, . . . ,dn such that di > pj for all Since 
each /i G I we have (/i, ...,/„) C I and V (/i, ...,/„) D V (I). The excess intersection of 
(/ii /2, ■ • • , fn) with respect to I is defined as the quasiprojective variety V (/i, . . . , /„) \V (I). 
We define the excess number i?, (I; c?i, . . . , d„) of an ideal X to be the number of isolated 
solutions in V (/!,...,/„) \V (I). 

In this paper, we consider the case when the polynomials /i,...,/„ are chosen to be T- 
generic forms (Definition [2| , which has applications in machine learning and ideal regression. 
When the generators of the ideal I C C[xo, . . . , x„] are all generic forms or all monomials, the 
excess intersection of (/i, . . . , /„) is a finite number of reduced points in P" whose cardinality is 
determined hy di, ... ,dn [D H] • The cardinality of this set of points is exactly the excess number 

(I; di,..., d„) =#V (/!,...,/„) \V(Z). 

At times it will be more convenient to work with the number Eo (I; di, . . . , d„) :— di ■ ■ ■ dn ~ 
Et {I;di, . . . ,dn) for this case. This number is classically known as the equivalence of an ideal 
[3] [Chapter 6]. We will see how (T; di, . . . , d„) and Eo {I;di, . . . ,dn) relate to the volume 
of a subdivided simplex. Our main contribution is a combinatorial proof of the theorem below, 
and algorithms using numerical algebraic geometry to compute excess numbers of any ideal. 

Theorem 1. Let X be an ideal ofC[xQ,..., Xn] generated by generic forms Bi, . . . , Bk of degree 
pi, . . . ,Pk. If fi, . . . , fn are X-generic forms of degree di, . . . , d„, then 

n—k 

E, (X;di, . . .,dn) +Pi ■ • -Pfc ^ (^{-if Vn-k-sVs^ = did2 ■ ■ ■ dn 

5=0 

where 'Dn-k-S is the sum of all square-free monomials in di, . . . ,dn of degree n — k — S and Pg 
is the sum of all monomials in pi, . . . ,pk of degree S. 
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The paper is structured as follows. We consider the case when Z is a monomial ideal, 
and show excess numbers equal mixed volumes of certain poly topes (Lemma |6]). By further 
restricting to the case when the ideal I defines a complete intersection that is also a linear space 
(though not necessarily reduced) , we do a mixed volume computation (Lemma |9]) to get an 
explicit formula for excess numbers. By Proposition [5] and Lemma |10| we immediately get the 
main result in Theorem [T] 

While also having applications to machine learning jH [2] , the motivation for this work came 
from Mike Stillman at the 2012 Institute for Mathematics and its Applications Participating 
Institution Summer Program for Graduate Students in Algebraic Geometry for Applications 
when he proposed the following problem. Given an ideal defining a curve, can we determine 
its equivalence? We give an answer to this question when the ideal is generated by generic 
forms and when the ideal is generated by monomials. Our computations were performed with 
the software bertini and Macaulay2 [H E]. In addition, we remark that our work has direct 
relations to computing Segre classes as seen in (6j [7l [8] . 

2 The Monomial Case 

The key idea to Theorem[l]is to cast our excess intersection problem in the language of combina- 
torial geometry and prove Proposition [5) Since the following proofs will use Newton polytopes, 
Minkowski sums, and genericity, we set up additional notation here. 

The Newton polytope of a form / will be denoted as Af (f). The standard n-simplex is the 
convex hull of the origin eoj a-nd the standard basis of unit vectors ei, . . . , e„ in M". The binary 
operation, Minkowski sum, will be denoted as "+". We take a generic form to be a homogeneous 
polynomial with full support and random coefficients. For example cix + C2y + c^z + c^w is a 
generic linear form in C [x, y, z, w] whenever the are chosen to be random numbers. For a 
more detailed definition of generic form, we refer to [T] . We now give a precise definition of a 
system of I-generic forms. 

Definition 2. Let I be an ideal of C [a;o, . . . , Xn] generated by forms Bi, . . . ,Bi whose respective 
degrees are pi, . . . ,pi. If the forms /i , . . . , /„ are defined by 

" /i 

/2 
fn 

then /i , . . . , /„ are said to be a system of I-generic forms with degrees di , . . . , (i„ whenever aij 
are generic forms with positive degree equal to di — deg Bj . 

With this definition, the degree of each fi is always greater than the degree of Bj for all i, j. It 
is important to note the difference between a form with generic coefficients, an I-generic form , 
and a generic form. A form with generic coefficients is a homogeneous polynomial with random 
coefficients that does not necessarily have full support. For example, whenever Ci, C2 are random 
coefficients, the form cix + C2y € C [x,y, z,w] has generic coefficients, but is not a generic form. 
An important fact is that an X-generic form need not have generic coefficients! For example, 
consider the ideal T = {x + y) of C[x,y, z,w\. Then, an I-generic form of degree 2 may be 
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written as 

{cix + C2y + C3Z + c^w) (x + y) = 
cix^ + C2y'^ + (ci + C2) xy + C3 {zx + zy) + C4 {wx + wy) . 

Since the coefficients of zw and zy are equal (as well as the coefficients of wx and wy), this 
I-generic form does not have generic coefficients. But, whenever I is a monomial ideal, an 
I-generic form will always have generic coefficients. This is key to proving Lemma |6] 

Remark 3. As mentioned in the introduction, a system /i,...,/„ of n X-generic forms of 
C [xq, . . . , Xn] has a variety consisting of V (I) and a union of finitely many points whenever 
the degrees of each fi are greater than the degrees of the generators of X. Also, if we take the 
T-generic forms /i, . . . , /„+; with T C C [xq, . . . , x^] and I > 0, then V (/i, . . . , /„+;) is equal to 
V(X) LI, 2J. 

Whenever I is a principal ideal, the excess numbers are easy to determine. 

Example 4. Let I ~ (B) be a principal ideal of the ring C [xq, . . . , Xn] with degi? = p. Then 
n I-generic forms may be written as 

fi = aiB, /2 = a2B, ...,/„ = a„B, 

with deg fi — di and being a generic form of degree di — p. To determine the excess number 
-E, (I; di, . . . , d„), we saturate (/i, . . . , /„) by I. Doing so, we conclude that the excess intersec- 
tion of (/i, . . . , /„) is defined by (ai, . . . , a„) and consists of finitely many points. By Bezout's 
theorem, it follows that (I; di, . . . , dn) = (di — p) (^2 — p) ■ ■ ■ {dn — p)- 

In general, determining the defining equations of an excess intersection is harder. But Lemma 
[6] allows us to determine excess numbers of monomial ideals directly using combinatorics. As 
a consequence, we can prove the following proposition, getting an explicit formula for excess 
numbers of certain ideals. 

Proposition 5. Let X he an ideal of C [xq, . . . , Xn] generated by generic forms x\^ , • • ■ , x^f!^ . If 
fi, . . . , fn are I-generic forms of degree di, . . . , d„, then 

E, {I;di, . . . ,d„) +Pi---Pk'Y^ (^{-if T^n-k-sVs^ = ^1^2 ■ ■ ■ dn 

s=o 

where V^-k-S the sum of all square-free monomials in di, . . . ,dn of degree n — k — 6 and Pg 
is the sum of all monomials in pi, . . . ,pk of degree 6. 

To prove Proposition [5] we first prove Lemma [6] This lemma shows that excess numbers of 
any monomial ideal can be computed using Bernstein's theorem |S] and mixed volumes. After 
this, we set up and do a mixed volume computation in Lemma [8] and Lemma |9] The results of 
this mixed volume immediately give us Proposition [5] 

Lemma 6. Let X C C [xo, . . . , Xn\ he a monomial ideal, and let fi, . . . , fn he X-generic forms 
of degrees di, . . . , (i„. Then the excess numher E, (I; di, . . . , d„) equals the mixed volume of the 
Newton polytopes Af (fi) ,■■ ■ ,M {fn)- 



3 



Proof. By Bernstein's theorem [2j the number of isolated solutions of /i, . . . , /„ in P" with non- 
zero coordinates equals the mixed volume of the Newton polytopes of /i , . . . , /„ whenever the 
forms /i, . . . , /„ have generic coefficients. Because I is a monomial ideal, the forms /i, . . . , /„ 
indeed have generic coefficients. So we need only show that an isolated solution of the forms 
fi, . . . , fn has nonzero coordinates. 

By way of contradiction, suppose y = (t/Oj- ■•!?/«) is an isolated solution of the system 
/i,...,/„ with zero for at least one coordinate. Without loss of generality we may assume 
Un = 0. Let ~ denote passing to the quotient ring C [xq, . . . , / (xn) = C [xq, . . . , 
Then (yo, yi, • ■ • , G P""^ is a solution to the system /i,...,/™ of X-generic forms of 

C [a;o, . . . , Xn~i]- But by Remark [s] there cannot be any isolated solutions to a system of n 
T-generic forms in the ring C [a;o, . . . , Xn-i], a contradiction. □ 

The key idea of this proof was to note the linear space defined by the generators of X 
are contained in the coordinate planes. Thus, these hyperplanes are "hidden" and Bernstein's 
theorem counts exactly the solutions we are interested in. So to determine excess numbers of 
monomial ideals, we determine mixed volumes. In general, mixed volume computations are 
complicated, but in certain cases there is hope for an explicit formula. The certain case we now 
consider is when I — (a;^^ , • ■ • , a:^*"' ) is generated by powers of unknowns. If /i , . . . , /„ are X- 
generic forms with degrees di, . . . ,dn, then we determine the excess number E, {I;di, . . . ,dn) by 
calculating the mixed volume of Af (/i) , . . . jN' {fn)- Recall that the mixed volume flTT [Chapter 
8.5] can be calculated by determining the coefficient of A1A2 • • • A„ in the polynomial defining 
the volume of the scaled Minkowski sum 

AiAA(/i) + --- + A„AA(/„). 

If we let & denote a certain simplex, then Lemma[8]shows that slicing (3 by a certain hyperplane, 
we can subdivide 6 into two convex polytopes ©o and ©i. We will choose the certain hyperplane 
so that ©1 = XiAf (/i) + • • • + AnA/" (/«). Since Vol 6 Vol ©0 + Vol ©i , we compute our desired 
mixed volume by determining the coefficients of Ai • • • A„ in Vol©o and Vol© (Lemma |9]). 
To elucidate our ideas we consider the following example. 

Example 7. Let X = {x^'-,X2^) be an ideal of the ring C [xq, Xi, 0:2: 2:3] and let /i,/2,/3 be 
X-generic forms of degrees di,d2,d^. Then the Newton polytope of fi will be the convex hull of 
two tetrahedra. Specifically, the Newton polytope of fi is the convex hull of the following eight 
points in M"^ of which six are vertices oi N {fi): 

(pi,0,0) (d„0,0) -pi,0) (pi,0,d, -pi) 

(0,p2,0) (0, dj,0) (P2,rfi-P2,0) (P2,0,di -P2) • 

To avoid confusion with points in P", we describe points in M"^ as (ui,U2,U3) rather than 
{xi, X2, X3). We may also describe the Newton polytope J\f {fi) by its 5 supporting hyperplanes: 

• the 3 coordinate hyperplanes, 

• the hyperplane defined by ui + -^2 + ^3 — rf^, and 

• the hyperplane defined by + ~ 1- 

The normal vectors of the 5 hyperplanes supporting M {fi) are the same for every i. Indeed 
they are the standard unit vectors ei, €2, £3, the vector ei + £2 + £3, and the vector + ^£2- 
By standard polytope theory [11 [Proposition 7.12], it follows that the scaled Minkowski sum 
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AiA/" (/i) + A2A/' (/2) + AsA/" (/a), has the same 5 normal vectors as those of A/" (/i). Indeed, the 
supporting hyperplanes are 

• the 3 coordinate hyperplanes, 

• the hyperplane defined by ui + U2 + W3 — (Aidi + A2(i2 + ^ad^)- E^i^d 

• the hyperplane defined by ^^ + ^^ — (A1 + A2). 

Now note that four of the five fiyperplanes are defining facets of a simplex whose volume is 
^ (Aidi + A2(i2 + Asda)'^. The fifth hyperplane subdivides the simplex as seen in Figure 1. By 
subtracting the volume of the white figure from the volume of the simplex, we attain the mixed 
volume by considering the coefficients of A1A2A3 in the difference. In this example, we would 
find the excess number satisfies E, (x^'^ a;2^; di, ^2, d^) +P1P2 {di + d2 + ds — pi — P2) — did2d^. 



Figure 1: 

We now precisely define the polytopes 6, 60, 61- Let D = Aidi + • • • + A„(i„ and A = 
Ai + • • • + A„. Then 6 is defined as the n-simplex whose ?i + 1 vertices are the origin and Dei. 
Moreover, the volume of © equals D"/nl which has a term did2 ■ ■ ■ dnM^2 • • • A„ when expanded 
out. If we define the hyperplane hp hy ^ + • • • + ^ — A, then hp slices 6 into two convex 
polytopes ©0 containing the origin, and ©1. We will prove that ©1 is the scaled Minkowski sum 
AiAA(/i) + ... + A„AA(/„). 

Lemma 8. Let /i,...,/„ be I-generic forms with degrees di,...,dn. If T — (x^^ , . . . , x'j^'' ) , 
then ©1 = AiAA (A) + • • • + A„AA (/„). 

Proof. We will show that the polytopes ©1 and XiAf (/i) + • • • + XnN' (/„) have the same n + 2 
supporting hyperplanes so that they must be equal. Specifically, we show that the supporting 
hyperplanes are 

• the n coordinate hyperplanes, 

• the hyperplane h^ defined by mi + • • • + u„ — Z?, and 

• the hyperplane hp defined by ^ + • • • + — A. 

The first n + 1 hyperplanes support the simplex ©. The last hyperplane hp slices © into two 
polytopes one containing the origin and a second which is by definition the polytope ©1. Now 
because XiAf (fi) has n + 2 supporting hyperplanes consisting of 

• the n coordinate hyperplanes, 

• the hyperplane defined by ui + • • • + u„ — A^d^, and 

• the hyperplane defined by + ' ' ' + ]J ~ -^i • 

By standard polytope theory [11 J [Proposition 7.12[, it follows XiAf (/i) + ■ ■ • + A„A^ (/„) has the 
desired supporting hyperplanes. □ 
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With Lemma [8] we are able to calculate the mixed volume by determining the coefficient of 
A1A2 • • • A„ in a certain integral as seen in Lemma [9] 

Lemma 9. With the previous notation, Vol (©0) ^ polynomial whose coefficient of X1X2 ■■ ■ A„ 
equals pi ■■■pu J2s=o (("!)'' '^n-k-s ■ Vs^ ■ 

Proof. By Lemma |8] it follows that the volume of ©0 equals 



Vol (60) = 



with the bounds of each integral inside the brackets with respect to d^^ being {xi = 0) — > 
(xi = D — Xi^i — • • • — xi) and A denotes the simplex in fc-dimensional space with k + 1 vertices 
of eo, Api •£!,... , Apfc • efc. 

By using the calculus fact [/•■■/ dj.^dx„_i ' ' ' ~ ^ {D — — ■ ■ ■ — xi)'' and the bi- 

nomial theorem, we have 

Vol (60) = j,J{D-Xk xiYd^ 

= ^ ELo ((-1)' Qd"-' J{xk + --- + X,)' d^) 

S=Q\\ ^1 [n-k-Sy. (5+k)\' i 

with r = n — k. It is known how to integrate a linear form raised to some power over the 
simplex. So to get the last equality, we use ^12j [Remark 9], that says J [xk + • ■ ■ Xif di^ — 

Now, note that Lm (Ai, . . . , A„) := m! (the monomials in Ai, . . . , A„of degree m) is congruent 
to A™ modulo A^, . . . , A^. Similarly, note that is congruent to Lm (Aidi, . . . , A„(i„) modulo 
A^ , . . . , A^j . So we have 

Vol (60) = pi • • ■PkJ2's=n (Aidi, . . • , A„d„) • Lk+s (Ai, . . . , A„)7'5^ 

= Xl ■ ■ ■ Xn ■ Pi ■ ■ ■ Pk J2's=0 ii^^ f 'Dn-k-sVs 



The last congruence is shown by an easy combinatorial argument. □ 

With Lemma|9] we have provided the explicit formula in Proposition [s] Since I' = (x^*^ , ■ • • , ) 
is a complete intersection, by perturbing the defining equations of I' by an e-amount, we create 
an ideal T which is a complete intersection defined by generic polynomials of C [xq, . . . 
Since we have perturbed the defining equations of T by an e-amount the isolated solutions in 
the excess of (/i, . . . , /„) with respect to I will also be perturbed by an order of e-amount. The 
point being that the number of isolated solutions is preserved and these two ideals have the 
same excess number giving us the following. 

Lemma 10. Let Bi, . . . , Bk he generic forms of C[xq, . . . , Xn] of degree pi, . . . ,pk. If li — 
{xP\. . .^xP"), then E, ...,dn)^E, {I;di, . . . ,d„)- 

As an immediate result of Lemma [TO] and Proposition [5] we have Theorem [T] 
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Remark 11. We remark that the number (I; di, . . . ,dn) depends only on the Newton poly- 
topes A/" (/i) ,A/'(/„). For example, consider the ideals Ii = [x^,y^),X2 = {x^ , y'^ , x'^ y , xy"^) , 
2^3 = {x^,y^j x^y^) in the ring C [w, x, y, z]. All three of these ideals have the same excess inter- 
section when every di is greater than 4 because the Newton polytopes of li-generic forms are 
the same for i — 1,2, 3. In particular, E, (1^; 5, 5, 5) = 44 for i — 1,2, 3. But if we consider the 
ideal = (^x^,y^^ + (xy), we find the Newton polytopes of J'-generic forms are different from 
those Newton polytopes of X^-generic forms. In particular, one can compute the excess number 
E, (J; 5, 5, 5) to be 65. 

3 Algorithms from Numerical Algebraic Geometry 

We have given a combinatorial description of excess numbers of monomial ideals in the first part 
of the paper and used this idea to give an explicit formula in Theorem [T] In the last part of this 
paper, we give algorithms that use homotopy continuation, an idea from numerical algebraic 
geometry to compute excess numbers of any ideal Z C C [xq, ■ ■ ■ , Xn\- 

We set the notation for this section. Let I C C [xq, . . . ,Xn] be an ideal generated by the 
forms Bi, . . . , Bi, and let T' G C [xo, ■ ■ ■ , Xn] be a monomial ideal generated by terms Aj where 
j is over some index set. Let /i, ...,/„ be I-generic forms, and let /(,..., be I'-generic forms 
such that di = deg/i = deg//. We implicitly assume that di is greater than deg Bj and deg Aj' 
for all possible 

Our algorithms will construct two homotopies, called h^pp and hjte, that take the isolated so- 
lutions of the Z'-generic forms /(,..., as start points and tracks them to solutions of /i, . . . , /„ 
giving bounds on E, (Z; di, . . . ,dn). In the first algorithm, the monomial ideal T' is constructed 
so that E, {I';di, . . . ,dn) > (X; di, . . . , By doing a numerical membership test |10| 
[Chapter 15], we will determine (X; di, . . . , c?„) and isolated solutions of V (/i, . . . , /„) \V (I) 
explicitly. In the second algorithm, the monomial ideal X' is constructed to give lower bounds 
of (X; di, . . . , dn) instead. But by iterating the second algorithm, we have a probabilistic way 
to make this bound sharp and compute all isolated solutions of V (/i, ... , /„) \V (X) explicitly. 
The hupp~homotopy gets its name because it produces an uppGv bound of E, (X; di, . . . , dn) prior 
to a membership test. The hitc-homotopy gets its name because several iierations can produce 
sharp lower bounds of E, {I;di, . . . , (i„) after a membership test. 

3.1 Algorithm one and the hupp-homotopy 

We now give a definition of the hupp-homotopy and prove that it does indeed provide an upper 
bound of Et (X; di, . . . ,dn) prior to a membership test. 

Definition 12. Let Bi,...Bi £ C[xo, . . . ,Xn] be forms such that Bj = J2k^3-h with Aj^k 
being a term. (^Note that Aj,k has both a coefficient and a monomial associated to it.) To ease 
notation, let Aj = [Aj^i, . . . , Aj^kj] be a row vector, Oi^j be a row vector of kj different generic 
forms, and be a row vector of a generic form repeated kj times. Define the hypp-homotopy 
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as hupp {t;di, . . . , dn) '■— 
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with the degrees of the generic forms of a^j and /3i.j chosen so that hupp (t,di, . . . ,d„) is a 
system of n forms of degrees di, . . . ,dn whenever t is a random complex number. We denote 
the start points of hupp(t;di, . . . as iSh^pp and take them to be the isolated solutions of 
hupp (1; di, . . . , c?„). Denote the end points of hupp(t; di, . . . ,dn) as Th^pp . 

With this definition, we have when t = 1 that hupp (i; di, . . . , dn) is a system of n I'-generic 
forms /{,..., of degrees di, . . . ,dn- On the other hand, when t = Owe have hupp (t; di , . . . , (i„) 
is a system of n I-generic forms of degrees di , . . . , d„ . By the fundamental theorem of parameter 
continuation of isolated roots [10 [Theorem 7.1.6[ it follows that Th^pp contains all isolated 
solutions of /i, . . . , /„. In particular, this proves Theorem 13 because #:S'hu — #^hu ■ 



Theorem 13. Let X C C[xq, . . . ,Xn] be generated by the forms Bi, . . . Bi such that Bj = 
'^if. Aj^k with Aj^k being a term. If we let X' be generated by Aj_k, then 

E,{I';di,...,d„)>E,{I;di,...,d„). 

Moreover, the parameter homotopy hupp (i; di, . . . , d„) has endpoints Ih^pp containing all isolated 
solutions of f[, . . . , . 

Now that we have the theorem, we present our algorithm. 

Input: Natural numbers di, . . . , d„ and generators Bi, . . . , Bi of an ideal I in C [xq, . . . , x„ 
such that Bj = Aj^k and Aj^k is a term. 
Output : The excess number E, {I;di, . . . , dn). 
Step 1 : Construct the the hupp-homotopy hupp (t; di, . . . , d„). 



Step 2: Solve the start system hupp (1; di, . 



[/(,..., f'n] and compute S'h„ 



Step 3: Use the hupp-homotopy to determine Thupp and an upper bound of E, (I; di, . . . , d„). 
Step J^: Use a numerical membership test to determine E, (I; di, . . . , d„) and the isolated solu- 
tions of hupp (0; di, . . . , d„) := /i, . . . , /„. 



For this algorithm, we assume in Step 2 that the excess intersection of a monomial ideal can 
be determined. While this may not always be true in practice, the first section has given some 
tools to address this case. We now give an example where I defines the twisted cubic. 

Example 14. Let the ideal Z C C [x, ?/, w\ be generated by the forms 

Bi ^ — yw, B2 = yz — xw, B^ = — xz, 
and suppose we want to calculate E, (I; 3,3,3). To run the first algorithm, we input di = d2 = 
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— 3 and Bi,B2, B3. In Step 1, we determine h^pp {t;di, . . . , dn) 



ail ai2 
0-21 0,22 
131 0,32 



ai6 

026 
036 



+ {l-t) 



hii hii bi2 bi2 bi3 bi3 

^21 ^21 &22 ^22 ^23 ^23 
^31 ^31 ^32 ^32 &33 ^33 



z 

—xw 
—xz 



The forms fly and hij are generic linear forms of ( 



[x, y, z, w]. Once we have solved the system 



hupp (1; 3, 3, 3) in Step 2, we path track in Step 3 to calculate Thupp giving an upper bound ^Th^ 
of E, (I; 3, 3, 3). In Step 4, we use a numerical membership test to conclude E, (I; 3, 3, 3) 
10. Indeed, if 
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then we find the ten excess points are si, S2, S3, S4, S5, sg, sy, ss = sq, sg = sg, sio = s?: 



Si S2 S3 S4 

a; -6.1999 -0.2081 -1.0024 -0.1530 
y 5.9766 0.5979 3.1208 0.3771 

z -2.3702 -2.1386 -5.1077 -0.6183 
w 1 1 1 1 



S5 

X -.6493 + 1.4057i 
y .4134- 1.4061i 
z -1.1267 + 0.3173i 
w 1 



0.4713 - 0.0461i 
0.2603 - 0.5271j 
-0.9278 + 0.1923i 
1 



S7 

2.9076 + 0.0384i 
-1.0341 + 1.7553i . 
-0.7082 - 1.2392i 
1 



3.2 Algorithm two and the hitc-homotopy 

While our first algorithm was a probability-one algorithm to determine the excess number 
E, {I;di, . . . , dn), the second algorithm is probabilistic. The algorithm uses the hjtc-homotopy 
to compute a lower bound of excess numbers. By iterating this algorithm, the lower bounds can 
become sharp. 

Definition 15. Let Bi,...Bi £ C [xq, . . . , x„] be forms that generate I and Ai,...,Ai be 
monomials that generate I' such that degi?j = deg Aj. The hjtc-homotopy is defined as 
hitc {t;di, . . . ,d„) := 



an ai2 ■ 


■ an 


/ 


' Ai - 




" Bi ' 








t 




+ 7(1-0 




] 




a,ni 


\ 


. . 




. Bi _ 
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with the degrees of the generic forms equal to deg/i — degAj. We denote the start points of 
hite as Site and take them to be the isolated solutions of hitc {l;di, . . . ,dn) and denote the end 
points of hitc {t; di, . . . , dn) as Titc- 

With this definition, we have when t = 1 that hjto (t; c?i, . . . , d„) is a system of n I'-generic 
forms /{,..., /4 of degrees di, . . . , (i„. On the other hand, when t = we have hjtc {t;di,..., dn) 
is a system of n I-generic forms of degrees di, . . . , c?„. While easy to set up, the hjte-homotopy 
is unable to use the fundamental theorem of parameter continuation of isolated roots [TO] [The- 
orem 7.1.6] so Thit„ does not necessarily contain all isolated solutions of /i, . . . , /„. However, 
after doing a membership test, we can determine some points in T^.^^ are isolated solutions 
of /i, . . . , /„. So what we have is a lower bound on {T;di, . . . ,dn). But, by iterating this 
honiotopy, we can find more isolated solutions and give a better lower bound. 



Input: Natural numbers c?i, . . . , dn, generators Bi, . . . , Bi of an ideal I C C [xq, . . . , Xn], mono- 
mials Ai,A2,...,Ai such that degAj = degi?j, and a (possibly empty) set W of isolated 
solutions of /i, . . . , /„. 

Output : A set Wuc containing W of isolated solutions of /i, . . . , /„, and ifWuc a lower bound 
for the excess number (I; c?i, . . . , rf„). 

Step 1 : Construct the hjto-homotopy hjto (t; c?i, . . . , d„) and track start solutions S'itc to target 
solutions Tite. 

Step 2: Use a membership test to determine which solutions of Tito are isolated and set Witc to 
be the union of W and isolated solutions of Tito. 

Step 3: Output Wuc and #Witc OR repeat steps 1 — 3 by making a different choice of 7 in the 
hito-homotopy. 



To do numeric calculation, we often carry out computations in a straight line program. In 
this regard, the hitc-homotopy would be effective because we need only making changes in the 
constant 7 to get new solutions. By taking different choices of 7 in the hjtc-homotopy we were 
able to produce the following example. 

Example 16. If we take Ai — z^^A^ = yz,A^ — , then we have the excess number 
E, (I'; 3, 3, 3) — 7. Next, we construct the hito-homotopy as hjto {t;di,..., dn) — 
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bl2 


bl3 


21 




623 


33 


bs.2 


633 





" z2 - 




yz 











z 



with aij the same as in Example 

■s'l J ^2 J '^3 1 *4 1 ^2 J ^3 1 ^4- 
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+ 7(1-0 

We find the 7 isolated solutions of hito (1; di, ^2, (^3) are 



yz - 

y' 



yw 
xw 
■ xz 



X 

y 

z 
w 



-8.4814 
8.2976 
-2.9043 
1 



-.0354-1- .78681 
.1201 -I- ~.7U6i 
-.9638-F .1650i 
1 



.8876 + .0702i 
.3006 - .5880i 
-1.2929 .2635i 
1 



-.3053 + .4774i 
.3779 - .7007i 
4.8276 .6092i 
1 



By talcing 7 to be different random complex numbers and keeping bij fixed, with 4 iterations, 
we were able to find that E, {X; di, . . . , dn) has a lower bound of 10. By the previous subsection, 
we know that this lower bound is actually sharp. 
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We comment that with certain choices of Ai,A2, A3 defining I, it can happen that E, (I'; 3, 3, 3) 
is greater than, equal to, or less than (I; 3, 3, 3). So one may be tracking too many paths, 
too few paths, or perhaps luckily the right number. Open questions remain about for which 
choice of monomials Ai^ . . . , Ai yield the best computational results. In addition, how should 
we choose 7 to guarantee new solutions will be found as we iterate; and how can we verify that 
our lower bound has become sharp are also interesting questions. These questions will remain 
for future work, and their answers may depend heavily on the context of the problem. 

Remark 17. We remark that the hjtc-homotopy need not have had the Aj be monomials. Any 
choice of a form Aj whose degree equals Bj could have been used. However, in this section, 
we have made the assumption that excess intersections of monomial ideals can be computed 
effectively, as we saw combinatorics can be used to understand excess numbers of monomial 
ideals. 

To conclude, we have shown that determining excess numbers of monomial ideals can be 
reduced to computing a mixed volume in certain cases. With this idea, we are able to provide 
an explicit formula for excess numbers of ideals with generic generators. We presented two 
algorithms using numerical algebraic geometry to determine excess numbers of any ideal. We also 
demonstrated that these algorithms have successfully lead to the calculation of excess numbers 
of an ideal defining the twisted cubic. We believe that ideals from combinatorics defined by 
sparse forms over many unknowns can have their excess numbers computed effectively using the 
hupp-homotopy, while the hite-homotopy may be used in any situation. 
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